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\ The purpose of this brief report is to present some results of our on-going project 

on the asymptotic behaviour of braneworld-typc solutions on approach to their pos- 
sible finite 'time' singularities. Cosmological singularities in such frameworks have 
served as means to attack the cosmological constant problem (see 1 and references 
therein). The main mathematical tool of our analysis is the method of asymptotic 
splittings introduced in Ref. 2 
I/"") ' Below we study a model consisting of a 3— brane configuration embedded in a 

five dimensional bulk space with a scalar field being minimally coupled to the bulk 
and conformally coupled to the fields restricted on the brane. The total action is 
• taken to be S to tai = S bu ik + S brane , where 
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= J d+xdY^E (Jjjs - | (V0) 2 ) , S b rane = -J xjgiftf) , at Y = n 



with Y denoting the fifth bulk dimension, k 2 = M t 3 , M* being the five dimensional 



o 

, Planck mass and f((j>) is the tension of the brane depending on the scalar field (p. 

We assume a bulk metric of the form ds 2 = a 2 (Y)ds 2 + dY 2 , where ds 2 is the four 
Q^- dimensional flat, de Sitter or anti-de Sitter metric. Then varying the above action 

| we obtain the field equations: 
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where k — 0, 1 or —1, and H" 1 is the de Sitter curvature radius. Assuming further 
that the unknowns are invariant under symmetry and solving the field 

equations on the brane we may express the solution in the form 

a'(n)--f/WK)) a (y;), ^(n) = //( ^ )} - (3) 

We now apply the method of asymptotic splittings to look for the possible asymp- 
totic behaviours of the general solution. Setting x = a, y = a' , z = , where the 
differentiation is considered with respect to T = Y — Y s (Y s being the position of 
the singularity), the field equations ((2]), become the following system of ordinary 
differential equations: 

x = y, y = —/3Az 2 x, z = —4yz/x, (4) 
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where A = k|/4. Hence, we have the vector field f = (y, —(3Az 2 x, — Ayz/x) J . Equa- 
tion {1} does not include any terms containing derivatives with respect to T; it 
is the constraint equation of the above system. In terms of the new variables, the 
constraint has the form 

y 2 /x 2 = A(3/iz 2 + kH 2 /x 2 . (5) 

Substituting the forms (x,y,z) — (aT p , r yT q , ST r ), with (p,q,r) 6 Q 3 and 
(a, 7, 5) € C 3 — {0}, in the dynamical system we seek to determine the possible 
dominant balances in the neighborhood of the singularity, that is pairs of the form 
B = {a, p}, where a = (a, 7, S) and p = (p, q, r). For our system we find: 

Bi = {(a, a/4, y/Z/Ay/Aff), (1/4, -3/4, -1)} (6) 
S 2 = {(a,a,0),(l,0,-1)} (7) 
B 3 = {(a, 0,0), (0,-1,-1)}. (8) 

Since @ is a weight-homogeneous system, the scale invariant solutions given by 
the above balances are exact solutions of the system. The balance B\ satisfies the 
constraint equation {2} only for k — 0, corresponding thus to a general solution for 
a flat brane, whereas B2 corresponds to a particular solution for a curved brane 
since it satisfies Eq. ([5]) for k ^ and a 2 = kH 2 . Finally the balance B3 represents 
a static universe conformal to Minkowski space and will not be analyzed further. 

Next we calculate the Kowalevskaya exponents, i.e., the eigenvalues of the matrix 
given by K, — D{(&) — diag(p); for B\ we find that spec(/C) = { — 1, 0, 3/2}, whereas 
for B2, spec(/C) = {—1,0,-3}. These exponents correspond to the indices of the 
series coefficients where arbitrary constants first appear. The —1 exponent signals 
the arbitrary position of the singularity, Y s . Since we have two non-negative integer 
eigenvalues the solution we are constructing will be a general solution (full number 
of arbitrary constants). 

Let us now focus on each of the two possible balances separately and build 
series expansions in the neighborhood of the singularity. For the first balance, we 
substitute in the system (J4j) the series expansions x = TP(a+ E^Cj-T 7 '/"), where 
x = (x,y,z), Cj — (cji,Cj2,Cj3), sis the least common multiple of the denominators 
of positive eigenvalues (here s — 2), and we arrive at the asymptotic solution 

, = ^4 32 T 7/4 +-, y = x', z = -^T- 1 --i^= C32 T^ + .... (9) 
7 iy/A lotyJAfi 

The last step is to check if, for each j satisfying j/s = p with p a positive eigenvalue 
corresponding to an eigenvector v of the K, matrix, the compatibility conditions 
hold, i.e. v T • Pj = 0, where are polynomials in c,-, . . . ,Cj-i given by JCcj — 
(j/s)cj = Pj. Here the corresponding relation j/2 = 3/2 is valid only for j = 3 
and the compatibility condition indeed holds. We therefore conclude that near the 
singularity at finite distance Y s from the brane, the asymptotic forms of the variables 
are a — >0,a'— >oo,0'— > 00. This is exactly the asymptotic behaviour of the solution 
found previously by Arkani-Hammed et al in Ref. 1 
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However, the previous behaviour is not the only possible one. The second balance 
has two distinct negative Kowalevskaya exponents and we therefore expect to find 
an infinite expansion of a particular solution around the presumed singularity at 
Y s . Expanding the variables in series with descending powers of T, in order to meet 
the two arbitrary constants occurring j = —1 and j = —3, and substituting back in 
the system ((4]) we find the forms 



Therefore as T — > 0, or equivalently as S = 1/Y — > oo, we have that a — > oo, 
a' — > oo and <£>' — > oo. 

We thus conclude that there exist two possible outcomes for these braneworld 
models, the dynamical behaviours of which strongly depend on the spatial geometry 
of the brane. For a flat brane the model experiences a finite distance singularity 
through which all the vacuum energy decays, whereas for a de Sitter or anti-de 
Sitter brane the singularity is now located at an infinite distance. We can choose 
the coupling such that to allow only for that behaviour met in flat solutions and, 
in fact, the particular form of the coupling used by Arkani-Hammed et al in 1 is the 
only choice to make this possible. This easily follows by using equations (J3]) and 
solving the Friedmann equation (flj on the brane for kH 2 , i.e. 



Clearly then k is identically zero if and only if f'(4>)/f(</)) — (2/3/3)k 5 , or equiv- 
alently, if and only if f(cf>) oc e^ 2 ' 3 / 3 * )K5 ^ (Arkani-Hammed et al in 1 have [3 = 3). 
By working with other couplings we can allow for non-flat, maximally symmetric 
solutions to exist and avoid in this way having the singularity at a finite distance 
away from the position of the brane. 
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